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ABSTRACT 

The galaxy cluster 1E0657-56 (z = 0.296) is remarkably well-suited for addressing outstanding is- 
sues in both galaxy evolution and fundamental physics. We present a reconstruction of the mass 
distribution from both strong and weak gravitational lensing data. Multi-color, high-resolution HST 
ACS images allow detection of many more arc candidates than were previously known, especially 
around the subcluster. Using the known redshift of one of the multiply imaged systems, we deter- 
mine the remaining source redshifts using the predictive power of the strong lens model. Combining 
this information with shape measurements of "weakly" lensed sources, we derive a high-resolution, 
absolutely-calibrated mass map, using no assumptions regarding the physical properties of the under- 
lying cluster potential. This map provides the best available quantification of the total mass of the 
central part of the cluster. We also confirm the result from Clowe et al. (2004, 2006a) that the total 
mass does not trace the baryonic mass. 

Subject headings: cosmology: dark matter - galaxies: clusters: general - gravitational lensing — 
galaxies:clusters:individual:lE0657-56 



1. INTRODUCTION 

The cluster of galaxies 1E0657— 56 is one of the hottest 
most X-r a.v luminous clusters known. Since its dis- 
covery by iTucker et al.l |1995|), it has b een the subject 
of intense and ongoing research LMehlert et al.1 120^ 
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partic ular, Chandra observations by IMarkevitch et all 
()2002tl revealed the cluster to be a supersonic merger 
in the plane of the sky with a textbook example of a 
bow shock, making this cluster a unique case in which to 
study hydrodynamical properties of interacting systems. 
The optical images show that the cluster has two dis- 
tinct components, and the X-ray analysis reveals that the 
lower mass sub-cluster has recently exited the core of the 
main cluster with a relative velocity of 4500^8^™ km s~^. 
Although this relative velocity appears to be unusually 
large, an analysis of cosmological simulations demon- 
strates that it is well within the pred icted range of the 
currently favored cosmological model l)Havashi fc Whiti^ 
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Due to its unique geometry and physical state, this 
cluster is the best k nown system in wh ich to test the dark 
matter hypothesis l|Clowe et al.ll2003) . The observed off- 
sets between the weak gravitational lensing mass peaks 
and the X-ray gas component give the most direct evi- 
dence for the presence o f dark matter yet available . Using 
the same observations, IMarkevitch et al.l (|2004D placed 
upper limits on the dark matter self-interaction cross sec- 
tion. 

The goal of this work is to obtain a high-resolution, ab- 
solutely calibrated mass map with no assumptions on the 
physical properties of the underlying cluster potential. 
For this purpose we use HST ACS data, incorporating 
the gravitational lensing information from both multiple 
image systems (strong lensing) and from distortions of 
background sources (weak lensing). The superb spatial 
resolution delivered by ACS both increases the number 
density of background sources that can be used for weak 
lensing and reveals new strong lensing candidates, espe- 
cially around the subcluster. The joint strong and weak 
lensing analysis thus benefits greatly from the ACS data, 
enabling us to increase the spatial resolution and signal 
strength of the mass map in the core of the cluster. 

Encouraged by the success of the combined strong 
and weak lensing reconstruction method developed in 
.Bradac et al. (2005b) (hereafter Paoer I) and the results 
from applying the ni ethod to cluster RX J1347.5— 1145 
ijBradac et al.ll200'5aj) . we proceed to apply this method 
to 1E0657— 56. Because the strong lensing data are richer 
in this case (RX J1347.5— 1145 has not been observed 
with either the ACS or WFPC2 cameras prior to this 
study), we improve the method as described in the text. 

We perform the reconstruction in the following se- 
quence. Using ACS (multi-color where available) HST 
images, we identify the multiply imaged systems. Hav- 
ing a spectrosc opic redshift for one lensed system (from 
IMehlert et ZII2OOI . we then use the predictive power 
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of simple strong lens modeling to estimate the redshifts 
of the other systems. Using positions and redshifts of 
the strongly lensed images and the shape measurements 
of the "weakly" lensing sources, we perform a combined 
strong and weak lensing mass reconstruction, thereby sig- 
nificantly improving the constraints on the mass and po- 
sitions of the main cluster and the colliding subcluster of 
1E0657— 56. In terms of the structure of this paper, in 
section|21we describe the optical images used in this anal- 
ysis, the basic image processing and the extraction of the 
strong f i)2.1|l and weak f ii2.2() gravitational lensing data. 
We then infer the mass distribution of cluster 1E0657— 56 
from these data in sectional following a demonstration of 
our methodology on suitable simulated data in section|31 
We discuss the possible sources of error in section O and 
summarize our conclusions in section |S1 

2. OBSERVATIONS AND DATA REDUCTION PROCESS 

ACS/WFC imaging of the cluster 1E0657-56 was car- 
ried out in Cycle 13 (proposal 10200, PI Jones) on 2004 
October 21 in two pointings with one and three different 
filters respectively. The two pointings are centered on 
the main cluster and the subcluster with a small overlap 
between them. The subcluster was observed in three dif- 
ferent filters (three orbits with F814W, one with F606W, 
and one with F435W), while the main cluster was ob- 
served only in the F606W filter (one orbit). 

The demands placed by the lensing analysis require 
special care w hen reducing the imag es. We use the 
Multidrizzle fKoeke moer et aLll2002 l) routine to align 
the images. To register the images with the astromet- 
ric accuracy needed for lensing analysis, we determine 
the offsets among the images by extracting high S/N 
objects in the individ ual, distortion corrected e xposures. 
We use Sextractor ("Berti n fc ArnoutsI [199^ and the 
IRAF routine geomap to identify the objects and calcu- 
late the residual shifts and rotation of individual expo- 
sures, which were then fed back into Multidrizzle. We 
use "square" as the final drizzling kernel and an output 
pixel scale of 0.03 arcsec; this is smaller than the origi- 
nal pixel scale of the ACS CCD to reduce the impact of 
resampling on the shape measurements. 

2.1. Strong lensing image identification 

Owing to its complex structure, the strong lens- 
ing analysis of the 1E0657— 56 cluster is complicated. 
Among the multiply- imaged systems, only the giant arc 
on the NW side of the main cluster's cD has a measured 
redshift (z = 3.24 lMehlert et alJl200j) . Corresponding 
images for multiply- imaged sources were identified by 
matching both morphologies and surface brightnesses in 
each of the available ACS bands (F435W, F606W, and 
F814W) and ground-based f ilters, BVR data f rom Mag- 
ellan and I-band from VLT ijClowe et al.ll2004|) . 

We confirm (based on photometry and morphology) 
that the 6 sy stems (labeled A-F i n Table ^ previously 
identified by iMehlert et alJ l)2001|) are indeed multiply 
imaged. We identify four additional systems (G-J) in 
the subcluster region, where none were previously known. 
All of these are identified in Fig. and image positions 
(corresponding to the peak surface brightness) are pre- 
sented in Table Unfortunately, for the greater part of 
the main cluster, only a single band ACS image is avail- 
able, which makes identification of additional candidates 



significantly more ambiguous. 

Using these identifications, we perform a parametrized 
strong lensing reconstruction. At this stage we are not 
generating a detailed strong lensing model, rather we 
use this parametrized model only to predict the red- 
shifts of the systems where spectroscopic redshifts are 
not available (all systems but A) and as the initial model 
guess for the subsequent strong and weak lensing recon- 
struction. For this reconstruction, we use only the im- 
age positions as constraints. The parametrized model 
consists of two non-singul ar isothermal ellipses (NIE) 
l)Keeton fc Kochaneklll99^ . for which the scaled surface 
mass density n is given by 

<0) = , ^ . (1) 

9 is defined counterclockwise from due west and h is re- 
lated to the line-of-sight velocity dispersion a through 
b oc a^. The components are centered on the southern 
cD galaxy of the main cluster and the brightest cluster 
galaxy (BCG) of the sub cluster respectively (denoted 
with white crosses in Fig. ^ it is not straightforward to 
identify the BCG of the main cluster, therefore we refer 
to the two main galaxies as the southern and northern 
cD). We allow the scaling b, core radii Vc, ellipticities |eg|, 
and position a ngles (j)^ of these to vary.^ Following the 
prescription of iKneib et al.l l)1996f) . we also include the 
20 brightest cluster members from the F606W-band (se- 
lected using color information from Magellan data) in the 
mass model. The galaxies are modeled as non-singular 
isothermal spheres with a line-of-sight velocity dispersion 
Cmcmb and core radius rcmcmb following 

Cmomb OC £ ^ , rc.niomb OC L ^ , (2) 

The proportionality constants are allowed to vary 
as well and below we quote the values of crmcmb 
and rcmcmb for a fiducial galaxy with a F606W- 
band magnitude mpeoew = 18. The best fit model 
for this system has values of {cr, rc, |eg| , </>g} = 
{1300 km s~\ 0.3', 0.2, 45°} for the main clus- 
ter, {cr,rc, |eg| ,0g} ^ {1000 km s~\ 0.3', 0.2, 280°} 
for the subcluster, and {o-incmb,i8, ?'c,moinb,i8} = 
{250 km s^^, 0.3'} for the cluster members. We also 
allow the redshifts of systems B-J to vary. The resulting 
best fit redshifts are given in Table To evaluate the 
angular diameter distances throughout the paper we 
assume the ACDM cosmology with ilni — 0.3, fl\ = 0.7, 
and Hubble constant Hq = 70 kms~^ Mpc~^. 

For objects close in projection to the main cluster, the 
redshift estimates are fairly reliable because the well- 
measured arc A is located there (we estimate an er- 
ror of Az = 0.2(1 -I- z)). However, the redshift es- 
timates for those close to the subcomponent are less 
well-determined, and are somewhat degenerate with the 
adopted mass of the subcluster. We postpone an ex- 
tension of the methodology described in the following 
sections that incorporates simultaneous solution for the 
unknown source redshifts to future work, but here note 

^ We define all ellip t icities e throughout this paper as in 
IBartelmann fc Schneidea l2n(!l^ with |e| = (l-/)/(H-/), / being 
the axes ratio. 
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TABLE 1 
The properties of the 
multiply-imaged systems used in 
this work. 





Ra 


Dec 


•^prcd 
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A 


104.62988 


-55.943798 








-0O.y4io44 




B 


104.63042 


-55.941474 


4.8 




104.63775 


-55.941851 




C 


104.63338 


-55.945324 


2.1 




104.64709 


-55.943575 




D 


104.63528 


-55.951836 


1.4 




104.640U8 


-55.95062U 




E 


104.64232 


-55.948784 


1.0 




104.65155 


-55.956671 




F 


104.64778 


-55.958284 


0.8 




104.56568 


-55.939832 




G 


104.56402 


-55.942113 


1.3 




104.56417 


-55.944131 






104.56293 


-55.939764 




H 


104.56133 


-55.942430 


1.9 




104.56189 


-55.947724 






104.56186 


-55.946114 




I 


104.56052 


-55.942930 


2.1 




104.56141 


-55.944264 






104.56909 


-55.946016 




J 


104.57025 


-55.944050 


1.7 



"The redshift of arc A was de- 
termined s p ectros copically by 
IMehlert et"aLl 120011^ . Other red- 
shifts are determined using the 
best-fit strong lensing model. 



that: 1) the parametrized model allows the generation of 
a set of source redshifts that, while model-dependent, are 
self-consistent; 2) the predicted redshifts are consistent 
with the colors of the objects. We adopt the redshifts 
estimated here for the combined strong and weak lensing 
reconstruction, and where possible estimate the system- 
atic error introduced as a result of this simplification. 

2.2. Weak lensing catalogs 

In the weak lensing analysis, we use the F435W, 
F606W, and F814W exposures for the subcluster, and 
F606W wherever multicolor ACS data are not available. 
We will only outline the main steps used for this anal- 
ysis, full details on how the weak lensing catalogs were 
generated will be given in a future paper dealing with the 
weak le nsing mass measure ments of 1E0657— 56 at large 
radius ijClowe et al.ll2006'^ in prep.). 

We correct galaxy images for the PSF anisotropy and 
P SF smearing c losely following the technique described 
in lClowe et al.l i2006bD . The proc edure is based on the 
KSB algorithm l|Kaiser et alJll995|) . in particular we use 
the modified IMCAT implementation. The KSB method 
is formally valid only in the weak lensing regime (i.e. at 
radii much larger than where multiply imaged systems 
form and where the distortions are small). Howe ver re- 
cent simulations from the STEP project ifHcvman s et al.l 
l2006j) show that this approach is also valid in the non- 
weak lensing regime (i.e. with shear values of I7I 
0.1 — 0.3) to the accuracy needed for a single cluster 
reconstruction. Lastly, in the vicinity of the highly elon- 
gated arcs we are dominated by the strong lensing signal. 



thereby minimizing the effects of signal dilution due to 
imperfect PSF correction. 

We select stars from the half-light-radius vs. magni- 
tude diagram and fit a fifth order polynomial to their 
measured cllipticities as a function of their position on 
the drizzled image. After the correction the rms of the 
stellar ellipticity components e i and £2 (as defined in 
iBartelmann fc Schneider! 1)200 If) ) changes from 0.017 to 
0.015 and from 0.004 to 0.003, respectively, and the mean 
is shifted to in both cases. The shear is then measured 
independently in all available filters and a weighted av- 
erage is calculated from those to arrive at our final, PSF 
corrected ellipticity estimates. We use the ~ 1' overlap 
between the two F606W pointings to test for systemat- 
ics in the PSF correction. The differences in the shears 
derived from the two exposures is consistent given the 
noise level. 

The ACS camera has a PSF that varies both spatially 
and temporally. While the approach described above ac- 
counts for the spatial dependence, it does not account 
for time variations. We investigate the PSF time de- 
pendence by measuring stellar shapes in individual ex- 
posures (four exposures were taken within each orbit). 
We conclude that the PSF is not stable (i.e. the global 
pattern of the PSF changes) even within a single orbit. 
However our detected signal is much stronger than the 
typical residuals due to imperfect PSF removal caused 
by the temporal variation; therefore we use the combined 
images to measure shapes. In addition, the PSF shape of 
an individual star changes as a function of its radius (see 
e.g. iHevmans ct al. 2005, Jec ct al. 2006); we therefore 
match the radius used to measure the stellar correction 
factors with that for the galaxies. Using simulated weak 
lensing data we identified a constant bias between the 
measured and estimated shear values; this bias was cal- 
culated for each image set separately and was introduced 
as a multiplicative shear calibration factor when produc- 
ing final catalogs. 

The color cuts applied to each catalog to remove clus- 
ter and foreground galaxies were color-color cuts esti- 
mated using photometric-redshift templates of galaxies 
at redshifts z <2. We used templates of all but extreme 
starburst galaxies, under the assumption that there are 
not likely to be many faint, low-redshift starburst dom- 
inated galaxies, particularly in a cluster. As such, we 
expect that unless there is an unknown population of 
dwarf galaxies in the cluster with colors much different 
from known populations, we will have removed all but a 
handful of outliers (extreme starbursts or dust obscured). 
We only use objects which were detected in more than a 
single band; in particular we reject all the objects that 
are detected only in the F606W pointing of the main 
cluster and were too faint for detection in the Magellan 
images. 

Unfortunately the observations at hand do not allow 
us to determine the photometric redshifts of the sources 
used for weak and strong lensing. We would need at 
least NIR data to reliably determine the redshifts, as 
most o f the sources are loc ated at z ^ 1. Therefore fol- 
lowing [Clo^^^Q l|2006b|) . we estimate the redshifts of 
the background sources by applying the color and magni- 
tude cuts that remove the likely foreground popula tion to 
the HD F-S photometric redshift catalog of Fontana et alJ 
l)1999l) . For the remaining sources, we average the ratios 
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Fig. 1. — The two F606W pointings of the main and sub component region of 1E0657— 56. Multiply imaged systems are marked and 
labeled (see also Table 0. White crosses denote the positions of the southern cD of the main and BCG of the subcluster. 



DdDds/ Ds (where I?d, and D^s are the angular di- 
ameter distances to the cluster, source and between the 
cluster and the source respectively) to generate a very 
crude estimate of the individual redshifts, but good av- 
erage redshifts per population (provided the redshift dis- 
tribution of the HDF-S is representative of that of our 
background sources). As a result all of the galaxies in 
a certain magnitude bin will be given the same redshift. 
The main reason of this approach, instead of calculating 
a mean for the whole catalog is the fact that the im- 
ages have different depths in different parts of the recon- 
structed area, and this method allows us to account for 
the sudden jumps in the mean background galaxy red- 
shift across the image boundaries, which would otherwise 
introduce spurious features in the reconstructions. 

This method is only approximate, therefore we test 
the reconstruction method using mock catalogs with only 



limited knowledge of the redshift of the source population 
in Sect. |31 The weak and strong lensing reconstruction 
performs well under these conditions. In addition, any 
significant dilution of the signal will be compensated by 
the strong lensing signal (if more than one strong lens- 
ing system is used). The cluster is also at a moderate 
redshift, therefore changing the sources from redshifts of 
z = ltoz = 1.5 would change the lensing strength by 
14%. In conclusion we expect the contributions due 
to errors in the assumed distribution of D^Dds/ Ds to 
be smaller than that expected from the galaxy ellipticity 
shot noise. 

Finally, we fit the non-singular isothermal sphere NIS 
model to the the averaged tangential ellipticities (et)- 
Both components generate large tangential shear signal. 
We center circular bins on the southern cD of the main 
cluster and the BCG of the subcluster. The resulting NIS 
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fit for the main cluster results in (Tnis — 1300 kms , 
a^cNis = 0'.4, and (Tnis = 1000 kms^^, Xc.nis = 0'.6 for 
the subcluster. This is an oversimplified model because 
the cluster is morphologically disrupted and we do not 
correct for the contamination of the main cluster profile 
in the sub cluster region (and vice versa). 

3. STRONG AND WEAK LENSING 
RECONSTRUCTION FROM HIGH RESOLUTION 
(ACS) DATA 

Our strong and weak lensing reconst ruction is per- 
formed using the method described in IPaoer J . The 
main idea behind the method is to describe the cluster 
mass distribution using a set of generic model parame- 
ters, which we choose to be the projected gravitational 
potential i/'fc measured on a regular square grid. The 
total is defined as 



(3) 



where Xe(V'fc) is the contribution from weak lensing, 
Xm(V'/c) from strong lensing and R{^k) is the regular- 
ization function with regularization para meter r?. We 
closely follow the methodology presented in lPaperl , with 
the few differences described below. * 

The most important difference is that we extend the 
formalism to more than a single set of multiple images. 
This may appear straightforward because we only need to 
sum up the contributions to Xu from each system. Given 
A^s multiply imaged systems having Ni images each, we 
define Xm ^-s 



Ns Ni 
,2 _ jTT 



= 2^ 2^ bm,^ 3 



(4) 



where S is the covariance matrix and bm,i = Om,i — 

Z{zs^i)d{9m,i) — /3s, i- Ps,i is the average position of 
the i-th source (calculated from the image positions 
and given the deflection angle a), and Z{zs_i) is the 
so-called "cosmological weight" function as defined in 
iBartel mann fc SchneideJ ll200l . 

In iPaper j we assumed that the covariance matrix is 
diagonal and indepe ndent of the lens m odel parameters. 
As discussed in e.g. IKochane'3 lj2004j) . this approach is 
not optimal and can lead to solutions that are biased to- 
ward models that predict high magnification at the im- 
age positions. With more sources at differing redshifts, 
this problem becomes more apparent because the method 
likely converges toward a solution with steep gradients 
in the surface mass density close to image positions (al- 
lowing for the "fine tuning" of the reconstruction). To 
overcome such unphysical solutions, we use a better ap- 
proximation for the error containing the magnification 
at the image position ^{Ora.i) and obtain the covariance 



matrix S 



M(^m,i) 



diag(as 1, (Tg 2)- Such a procedure 

is approximately corre ct, provided we are close enough 
to the true model fsee lKochanekll2?)0^ . The covariance 
matrix is still diagonal; however we argue that assuming 
it to be such in the source plane is in fact a better ap- 
proximation than assuming the covariance matrix to be 



® In figure 2 in IPaper H the factor of ^ 5 in finite differencing 
formula for k is wrong. All the calculations in all three papers do 



diagonal in the image plane, as sources are on average 
more circular than their lensed images. In addition, in 
practice multiple image constraints are satisfied nearly 
perfectly and exact values of errors on image positions 
are of lesser importance. 

We adapt the method to accommodate data regions 
of arbitrary size and shape, but the potential is recon- 
structed only in grid cells containing weak and/or strong 
lensing data. We tile the observed field with a regu- 
lar grid and determine which cells contain either weak 
lensing galaxies or multiple images. This can result in 
some cells being sparsely populated (holes in the data, 
edge effects), but the regularization ensures that poten- 
tial is properly reconstructed. Once the data-grid is de- 
termined, we ensure that each of the data grid points is 
surrounded by a 4 x 4 grid of points that are included 
in the reconstruction to calculate the scaled surface mass 
density k, shear 7, and the deflection angle a, using finite 
differencing, as described in Paper I. 

To find the x^"function minimum, we search for a so- 
lution of the following system of equations ^'g^^'''' — 0- 
This is in general a non-linear set of equations, which we 
solve in an iterative manner. We again linearize this sys- 
tem as described in iPaperT . calculating all of the non- 
linear terms (for the strong lensing term in particular 

these are now the ones containing /3s, i and fJ-{Om,i)) using 
the information from the previous iteration. We reg- 
ularize the solution using a "moving prior" , gradually 
updating our knowledge of the model by increasing its 
complexity and computing its likelihood relative to the 
previous (simpler) model. We begin with an initial model 

(k^"^ , 7^°^ , 72°^) computed on a relatively coarse grid, and 
then gradually increase the number of grid points, com- 
paring the resulting k map with that from the previous 
iteration, linearly interpolated onto the finer grid. In ad- 
dition to K we now also use the components of the shear 
71 and 72 in the regulariz ing sche me, to avoid the nu- 
merical effects described in lPaoei^ even more efficiently. 
The regularization constant was set after experimenta- 
tion with simulated data. When using many multiply 
imaged systems, the final reconstruction is even less sen- 
sitive to thes e choices than when only a single system 
is used (as in iPaoerT) : we discuss this finding in more 
detail in Sect. ]E\ 

With many multiply imaged systems, we must choose a 
sufficiently fine grid to be able to calculate the difference 
in the deflection angle between images of different sys- 
tems (if images from more than one system lie within the 
same grid cell, it is difficult to find a good solution). Ide- 
ally one would like to use adaptive grids for this purpose 
(having a relatively coarse grid where only weak lensing 
information is available and increasing the resolution in 
the vicinity of the multiple images). Unfortunately it is 
computationally difficult to track such a problem and is 
therefore beyond the scope of this paper, but will be a 
subject of future work. 

We tested all the improvements described above us- 
ing a high-resolution N-bod y simulation of a galaxy clus- 
ter byiSpringel et al.l l|2001fl . The catalogs are generated 
as inLPaDer I, although adapted to the higher-resolution 
data that we use here. The weak lensing simulated data 
are obtained by placing 1800 galaxies (giving a den- 
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Fig. 2. — Average tangential ellipticity (et) vs. projected radius (6) in radial bins centered on the southern cD of the main cluster {left 
pannel), and the BCG of the subcluster (right pannel) containing 100 galaxies per bin. The errors are obtained by randomizing the phases 
of the measured ellipticities, while preserving their absolute values. 100 realizations are performed. The dashed line shows the best-fit 
non-singular isothermal sphere (NIS) profile to the binned data (see text). It is plotted here for the average source rcdshift of weak lensing 
sources (2:) = 1.0. 



sity of 120 arcmin"^) randomly on a 3.8 x 3.8 arcmin^ 
field. The intrinsic ellipticities eg are drawn from a Gaus- 
sian distribution, with each component characterized by 
a = (Tes = 0.2. We then draw the redshifts of the back- 
ground sources, following .Brainerd et al. (1996), from a 
gamma distribution Pz{z) = z^/(2zq) exp(— z/zq) with 
zq = 1/3, giving a mean redshift of (z) = 3zo = 1, and 
a mode of Zmode = 2zo = 2/3. We reject the foreground 
objects, and use the background object redshifts when 
predicting the lensed shapes. However, when using the 
catalogs for the reconstruction, we assume no knowledge 
of the individual redshifts: instead all sources are as- 
sumed to have z = 1, in agreement with the poor esti- 
mates available for the actual data. The multiple -image 
systems are generated as described in iPaper 1 Four 
quadruply- imaged systems are generated at source red- 
shifts of Zs = {0.8, 1.0, 3.0, 5.0}. The initial model for the 

regularization is set to k^'^'' — ^[^^ = 72"'' = 0. 

The resulting reconstruction is given in Fig. |21 It is 
very encouraging to observe how well the overall shape 
(the agreement in the ellipticity and position angle of the 
main component, and the overall profile) and its mass are 
reconstructed. These successes would not be possible 
with weak lensing data alone, owing to the mass-sheet 
degeneracy (i.e. the reconstruction would suS^er from the 
degeneracy of the form k ^ k' = = Ak -|- (1 — A), whe re A 
is an arbitrary constant, see e.g. 'Bradac et al. 2004). In 
addition with the ACS-quality data, we arc starting to 
resolve the small-scale substructure of the potential. 

4. CLUSTER MASS RECONSTRUCTION OF 1E0657-56 



We now apply the above mentioned reconstruc- 
tion method to the strong and weak lensing data of 
1E0657-56. As described in Sect. O we use 10 dis- 
tinct multiply imaged systems and A'gai ^ 1900 "weakly" 
lensed galaxies (120 arcmin^^). We start with a 30 x 
30 pix^ grid for a 7 x 7 arcmin^ field oriented with respect 
to the negative RA coordinate (not all of the grid cells 
contain data and those that do not are excluded from 
the reconstruction) . We gradually increase the number of 
grid points (for reasons described in Sect.|3J) in steps of 1, 
with the final reconstruction performed on a 60 x 60 pix^ 
grid. In our initial model, we use a two-component NIE 
model with parameters as given in Sect. 12.11 to set the 
regularization. The possible dependence of the recon- 
struction on initial conditions is discussed in the next 
section. 

The resulting reconstruction is shown in Fig. ^ to- 
gether with X-ray surface brightness contours from the 
500 ks Chandra observation. The two cluster compo- 
nents are clearly detected. In addition, with superb 
strong lensing data (due to exquisite ACS resolution and 
multi-color information for the subcluster) we are start- 
ing to resolve multiple components within each individ- 
ual system. 

From this reconstruction, we measure the enclosed 
(projected within a cylinder) mass for this system. The 
main cluster has Mmain(< 250 kpc) = (2.8 ± 0.2) x 
IO'^^Mq, while the sub cluster has Msuh{< 250 kpc) = 
(2.3 ± 0.2) X IO^'^Mq (1 arcmin = 260 kpc at the cluster 
redshift) . The mass errors represent combined errors ob- 
tained by bootstrap resampling of the weak lensing cata- 
logs, performing a similar experiment for the strong lens- 
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(arcmin) 9. (arcmin) 



Fig. 3. — The reconstructed (left) and the original (right) surface mass density k of the simulated cluster described in section|21 The 
symbols in the left panel denote the image positions of the four multiple-image systems at Zg = {0.8,1.0,3.0,5.0} which we use for the 
reconstruction. The contour levels, which are the same in both panels, are linearly spaced with Ak = 0.38, starting at k = 0.4 for a fiducial 
source at infinite redshift, Zs — ► oo. 



ing data, and estimating systematic uncertainties from 
the mass differences obtained using different initial mod- 
els (see discussion in the next section) . We plot the inte- 
grated, azimuthally-averaged mass profiles M(< R) for 
the two components in Fig. We fit power-law mod- 
els {M{< R) oc i?"') to these profiles, and find loga- 
rithmic slopes of rii^main = 0.8 for the main cluster and 
"■i,sub — 1-1 for the subcluster. In terms of the surface 
mass density K(r) cx r"", this translates to rimain = 1-2, 
and risub = 0.9. Adopting an isothermal profile {n = 1), 
we find the line of sight velocity dispersions for the main 
and the sub cluster are cr,nain — (1400 ± 100) kms^^ and 
fsub = (1200 ± 100) kms^^ res pect ively; in good agree- 
ment with the analysis in Sect. However, as already 
mentioned in that section, given the size of the errors, the 
asphericity of the system, and the complicated geometry 
(the relative contributions of the main and sub cluster 
are difficult to disentangle), the velocity dispersion and 
slope measurements should only be used as guidelines. 

Because of this cluster's violent history, we can not ex- 
pect the assumption of isothermality to hold. Therefore 
(and for the reasons stated above) it is also non-trivial 
to compare o ur velocity dispersion measurements with 
the ones from iBarrena et aLl l|2002D . Their estimate for 
the main cluster agrees with ours, while the velocity dis- 
persion of a subcluster is much lower in their case. An 
iso lated lens with the lo w velocity dispersion measured 
by IBarrena et al.l l)2002D for the subcluster would have 
multiply imaged systems only at radii < 1. 3, a factor 
of ^ 8 below that observed. The influence of the main 
cluster can increase the radial range over which multi- 
ple images exist, but not to the extent observed here. 
We conclude that the true velocity dispersion of the sub- 
cluster is likely to be smaller than that derived from our 
mass model (due to the inclusion of mass from the main 
cluster), but larger than that presented bv lBarrena et al.l 



(200^. 

One of the important results of the work of lClowe"et 
|2004, 2006a) was the measurement of the offset between 
dark matter and baryons. We confirm these results: our 
strong and weak lensing offsets for the mass peaks (by 
measuring the peaks in ^-distribution) and the corre- 
sponding errors are given in Table |21 To estimate the 
significance of this offset we also measured the positions 
of the two peaks (for the main and the subcluster) in each 
of the 1000 reconstructions using the bootstrap resam- 
pled weak lensing catalogs (described below). In none 
of these realizations are the mass peaks coincident with 
the corresponding X-ray peaks. The 1-a error bars are a 
factor of 10 and 6 smaller than the offsets for the main 
and subcluster respectively. If these errors are accurate 
and Gaussian distributed, this result translates into 10-ct 
and 6-a discrepancies between the corresponding gas and 
the total mass peaks. These offsets directly demonstrate 
that dark matter is present and is the dominant mass 
component of this cluster. 

Lastly, we compare the projected gas mass (the details 
of modeling will be given in Randall et al. 2006) vs. the 
total mass within elliptical regions centered on the two 
peaks of the X-ray emission and the southern cD of the 
main and BCG of the subcluster, respectively (see Ta- 
ble The ratio of gas mass to total mass is smaller in 
the regions centered on the massive galaxies rather than 
the X-ray peaks, indicating that the gas has been (to 
some extent) stripped, and is no longer coincident with 
the total mass. In addition, we detect an extension in the 
total mass map to the NW of arc A. This detection can 
be partly attributed to the hot X-ray gas (see Fig.Q. 

5. POSSIBLE SYSTEMATIC EFFECTS 

To assess the reliability of the reconstructed map, we 
generate lOOO bootstrap resampled weak lensing catalogs 
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TABLE 2 

Offsets (with l-o- errorbars) of the mass peaks 

AND gas peaks FROM THE CORRESPONDING SOUTHERN 
CD OF THE MAIN AND BCG OF THE SUBCLUSTER 
(marked WITH CROSSES IN FlG.0. 





Ax([arcmin]) 


Aj/([arcmin]) 


Main component 


0.26 ±0.08 


-0.01 ±0.09 


Main component (gas) 


1.06 


0.84 


Sub component 


-0.05 ±0.10 


0.00 ±0.07 


Sub component (gas) 


-0.73 


0.10 



and perform reconstructions on each of these. To further 
test the rehabihty of the strong lensing data we create 10 
different reconstructions each time removing one of the 
multiply imaged systems we use. The resulting K-maps 
do not change substantially, the main features (i.e. the 
ellipticity of the main and the subcluster) seen in Fig. 0] 
remain in all of the reconstructions. 

We also perform a reconstruction using the weak lens- 
ing data only (bottom- left panel of Fig. |HJ|. From only 
weak lensing, we clearly detect both the main and sub 
clusters, but are unable to resolve individual substruc- 
tures within each with high significance. This failure can 
in part be attributed to the limited field size used here in 
contrast to Clowe ct al. (2006a) (our method at present 
does not allow us to use larger fields for the reconstruc- 
tion, adaptive grid methods need to be employed then). 
In addition, the mass estimates from weak lensing recon- 
struction are unreliable (due to the the mass-sheet de- 
generacy) ; the resulting potential does not reproduce any 
of the multiply imaged systems observed. Traditionally, 
to overcome the problem of the mass-sheet degeneracy 
and measure the mass from the weak lensing mass recon- 
struction, one assumes a vanishing surface mass density 
at the edge of the field or adopts parametrized model fit- 
ting. The former is not applicable here because the field 
used is too small. The latter will yield reliable mass esti- 
mates if the assumed profile is correct, but it is difficult to 
guess the appropriate model for such a morphologically 
disturbed cluster. Combining strong and weak lensing 
helps us to break this degeneracy. However, the lack of 
counter arc candidates and redshifts of multiply imaged 
systems limits our ability to provide even stronger con- 
straints on the mass distribution. 

The errors on mass estimates quoted throughout the 
paper are obtained by summing in quadrature the er- 
rors obtained from the reconstructions using bootstrap- 
ping of the weak lensing catalogs and removing individual 
strong lens systems. In addition, we estimate and add 
systematic uncertainties resulting from differences in the 
reconstruction when using different initial models (as ex- 
plained below). In principle, we should incorporate all of 
these effects simultaneously in the reconstructions, how- 
ever this is too computational intensive and beyond the 
scope of this analysis. 

5.1. Initial conditions 

Methods involving strong lensing mass reconstruction 
with many multiple imaged systems are subject to de- 
generacies in parameter space. It is difficult to search 
the parameter space because the function has many 



local extrema. In our case, there are pronounced, nar- 
row secondary minima where individual systems are re- 
constructed well by the model (and the others less so) 
and high "walls" in space where images happen to 
form close to the corresponding critical curves. Because 
weak lensing data are noisy, additional, though less pro- 
nounced, minima in created. Given that our re- 
construction method does not guarantee that we find a 
global minimum, in principle wc would need to check the 
whole parameter space. This search is unfeasible with 
any available algorithms, given that the number of pa- 
rameters (i.e., pixels in the mass map) is > 1000. 

Instead, we opt to investigate the use of many initial 
models. We try both extremes; i.e. models that are 
clearly in disagreement with the data (zero initial con- 
ditions, single mass peak centered on the position of the 
gas) and models that use alternative methods to recon- 
struct the data (two-NIE model, two-NIE plus cluster 
members model from strong lensing mode ling, weak lens- 
ing mass model from lClowe et al]l200 6a'). Although the 
details of the reconstruction depend on the initial condi- 
tions (as can be seen from comparing all but bottom-right 
K-maps in Fig.O, the main features (SE-NW elongation 
of the main cluster, two mass peaks close to the main 
cluster cDs, an elongation in the mass map toward the 
gas peak, E-W elongation of the subcluster) are indepen- 
dent of the particular initial model we use. 

The resulting contribution of the weak lensing data 
remains the same for all sets of initial conditions 
(X? =i 7000 for 2iVgai = 3800 data points). The strong 
lensing contribution is lowest when the two component 
NIE model is used (xm = 20) and highest when the zero 
model initial conditions (xm = 1000) are used. Although 
the absolute values of the x^ have little meaning (i.e. 
one can change the regularization parameter to obtain 
a lower value for x? because the weak lensing data is 
noisy), the relative values tell us which resulting recon- 
struction is indeed a better fit to the data (as the same 
regularization parameters are used in all cases and it is 
the fit to the high signal to noise strong lensing data that 
improves). 

We estimate the systematic uncertainty on our mass 
estimates by calculating the maximum difference in the 
resulting mass estimates when using different initial mod- 
els. We exclude the zero model in this calculation as it 
does not predict the strong lensing data reliably. The 
resulting K-map (see bottom-left panel in Fig. IBJ for this 
model is suppressed where only weak lensing data are 
available (due to the mass-sheet degeneracy) and pro- 
duce steep gradients close to arc A and a ^ 10% lower 
overall mass. However, other features in the reconstruc- 
tion (positions and shapes of the mass peaks for the main 
and the sub-cluster) remain unchanged. In all cases, the 
resulting mass peaks "move" away from the initial guess; 
showing that the regularization is sufficiently fiexible to 
reproduce the features in the data. Therefore, we are 
confident that our results are not biased as a result of 
the particular choices of regularization and are indepen- 
dent of the initial conditions to the extent presented in 
Fig. El 

5.2. Strong-lens system identification and redshifts 

Finally, the mass measurements presented above de- 
pend on the reliability of the redshifts and identification 
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used for strong lensing reconstruction. These (except for 
system A) a re no t secure and depend upon the model 
used in Sect. 12. l1 In particular, the redshifts of the sys- 
tems around the subcomponents (and consequently the 
measurements of its mass) can change because the sub- 
component mass reconstruction does not depend strongly 
on the modeling around the main component (where the 
arc with the known redshift is located). 

The parametrized model used in Sect. I^TI to estimate 
the redshifts of the strong lensing systems likely deviates 
from the truth model. We tested the dependence of the 
assumed redshifts by performing 100 reconstructions in 
which we modified all of the strong lensing system red- 
shifts except for system A. For each reconstruction these 
uncertain redshifts were changed by assuming a Gaus- 
sian scatter with tJz = 0.5 and requiring that all sources 
lie behind the cluster. The resulting mass and position 
errors were smaller (by a factor of ~2) than the ones 
quoted throughout the paper, and the mass and position 
measurements were unbiased (as expected from such an 
analysis). The mass measurements would change, how- 
ever, if we systematically increase (decrease) the redshifts 
of all the systems. For example, if we decrease the red- 
shifts of the subcluster systems such that the lowest red- 
shift system of the subcluster (G) is at Zg = 0.5 (and 
scale the others accordingly), then the mass increases by 
a factor of ~ Z{1.3)/Z{0.5) = 2. If we instead increase 
all the redshifts, pushing the highest redshift system (I) 
to z = oo), then the mass of the subcluster will decrease 
by l.G/Z(2.1) = 1.2. Note however that these are just es- 
timates, since the redshift estimates of the weak lensing 
sources also plays a role, and the changes will therefore 
be even smaller. 

As a further test, we also perform 10 reconstructions 
where we remove one strong lens system at a time from 
the analysis. In Fig. [3 we show the two that differ the 
most from the reconstruction that uses all the strong 
lensing systems. When we remove system I from the re- 
construction (right panel in Fig. [7)) the double peaked 
structure of the subcluster disappears, while the east- 
west elongation remains. For the main cluster we see the 
largest deviations from the original map when we remove 
system D from the reconstruction (left panel in Fig. Ej) . 
However, these results do not mean that these structures 
are not real, rather they illustrate some possible system- 
atics when dealing with sparse strong lensing data. We 
are confident that system D is indeed multiply imaged, 
and images 12-13, at least, are also probably multiply 
imaged. 

These tests demonstrate that even though the strong 
lensing system identification and their redshift determi- 
nations are important to the details of the reconstruction, 
the errors arising from misidentifications/redshift mises- 
timation are within those quoted in the previous section. 

5.3. Predictive power of the model 

The predicted redshifts and strong lensing system iden- 
tification can be tested using spectroscopy. We are plan- 
ning just such a survey and will therefore be able to test 
the fidelity of our model in the future. In addition, using 
Spitzcr IRAC data we identify a candidate high redshift 
object. While bright in all IRAC bands, the faintness 
of this galaxy in optical and near infrared bands indi- 
cates that it could lie at z > 6. It has two components 



and our predicted critical curve for this redshift (indis- 
tinguishable from the z — s- cx) critical curve) is exactly 
where it should be if indeed these sources are merging 
images of the same source; they are denoted with crosses 
and plotted together with the critical curve in Fig. |H1 
Furthermore, we identify a candidate counter image with 
similar colors in the IRAC data just outside the critical 
curve. All three need to be further analyzed and con- 
firmed (using spectroscopy and/or NIR data) before we 
include them in the modeling. 

We also identify possible counter-images for systems 
B, J, H, and I as predicted from our model, they are 
indicated in Fig. O If these and high-redshift candidate 
are confirmed, the predictive power of our model would 
further indicate that the strong lensing data we use is 
reliable enough for this work. 

6. CONCLUSIONS AND FUTURE WORK 

Massive and interacting clusters, while quite rare, are 
remarkably well-suited to addressing outstanding issues 
in both galaxy evolution and fundamental physics. How- 
ever, in order to study the mass distribution, methods 
relying on hydrostatic (X-rays) or dynamical equilibrium 
are ill-suited for such systems. 

We have applied a mass reconstruction method based 
on strong and weak gravitational lensing to the cluster 
1E0657— 56. We use deep, high resolution optical data 
to identify objects belonging to the same multiple-image 
systems. The same data are used to obtain weak lens- 
ing catalogs allowing us to obtain a strong-|-weak lensing 
mass map of the cluster core. Our main conclusions are 
the following: 

1. Using the combined strong and weak lensing mass 
reconstruction we derive a high-resolution, abso- 
lutely calibrated mass map; we get projected, en- 
closed mass Minain(< 250 kpc) = (2.8 ± 0.2) x 
IO^^Mq around the main and Msub(< 250 kpc) = 
(2.3 ± 0.2) X IO^Mq around the subcluster. 

2. We detect the main cluster peak and a distinct 
mass concentration at the subcluster position, both 
clearly offset from the location of the X-ray gas in 
the system (at 10-cr and 6-(t significance for the 
main and the subcluster respectively - see Fig. 0] 
and Tab. 121). 

3. The majority of the mass is spatially coincident 
with the galaxies, which implies that the cluster 
mass must be dominated by a relatively collision- 
less form of dark matter. The high resolution data 
allow us to significantly detect the shapes of both 
the main mass component and the subcluster with 
no prior assumptions on their positions or profiles. 

We show that the cluster 1E0657— 56 is a very efficient 
strong lens; the area enclosed within the critical curve for 
highest redshift s ources is comparable to that of A1689 
IBroadhurst et al.l ((^05) . Consequently, we plan to con- 
duct a pencil beam survey for high-redshift (z ^ 7) ob- 
jects. At present this is not possible, because we lack 
both reliable redshift information for the 1E0657— 56 
multiple image systems, and multi-color high-resolution 
imaging data in the main cluster region. The latter is 
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Fig. 4. — The F435W-F606W-F814 color composite of the cluster 1E0657— 56. Overlaid in red contours is the surface mass density k 
from the combined weak and strong lensing mass reconstruction (for the purpose of this plot we recalculate the final K-map from top-left 
panel of Fig.|n]on a finer grid). The contour levels are linearly spaced with Ak = 0.1, starti ng at k = 0.5, for a fiduc ial source at a redshift 
of Zs oo. The X-ray brightness contours from the 500 ks Chandra ACIS-I observations IMarkevitch et ani200g) are overlaid in white. 
North is up and East is left, the field is 4. 9' X 3.2', which corre sponds to 1300 X 830 kpc'^ at the redshift of the cluster. The color composite 
was created following the algorithm from lLupton et^ H200I) . 



TABLE 3 

Comparison between the gas mass and the total mass for four elliptical regions 
centered on the maxima of the gas density map and the southern cd / bcg of the 

main / SUBCLUSTER. Ra AND Rf, INDICATE SEMI-MAJOR AND SEMI-MINOR AXES [Ra IS ORIENTED 
W.R.T. RA) OF THE REGION CONSIDERED. LAST ROW GIVES THE CORRESPONDING NUMBERS FOR 
THE ENTIRE ACS FIELD (RECTANGLE). 





-Ra X [kpc] ( [arcmin] ) 


A/gas [10^3 Af0] A/total 


[1013 Af0] 


Afgas /A/total 


Main (gas peak) 
Main (S cD) 


125 (0.47) X 250 (0.59) 


2.0 ± 0.2 
1.4 ± 0.1 


10.8± 0.6 
16.3 ± 0.9 


0.19± 0.03 
0.09 ± 0.01 


yub (gas peak) 
Sub (BCG) 


80 (0.31) X 80 (0.31) 


0.42 ± 0.04 
0.18 ± 0.02 


2.1± 0.2 
4.3 ± 0.2 


0.20 ± 0.06 
0.04 ± 0.01 


ACS field 


1300 (4.9) X 830 (3.2) 


13 ± 1 


91 ± 10 


0.14 ± 0.03 



needed to identify counter images for the strong lens sys- 
tems, and so improve the rehabihty of the mass map to 
the south-east of the main cluster. Once these data are 
obtained (HST Cycle 15), we will improve our method- 



ology to include an adaptive grid approach, which will 
greatly improve the constraints on the location of criti- 
cal curves for these high-redshift sources. 
The prospects for such a search are however very 



The mass distribution of 1E0657— 56 



11 



o 0.4 




Fig. 5. — The integrated mass profile of 1E0657— 56. The profile 
was determined by measuring the enclosed mass in cylinders, cen- 
tered on the southern cD of the main cluster (dots), and the BCG 
of the subcluster (crosses). We fit power law profiles to the data 
M(< R) oc R"' , the results are shown for the main (solid line) and 
for the subcluster (dashed line). 



promising, as we already have two strong candidate high 
redshift objects. One is the red arc B, which is con- 
strained to he at 3.24 < z < 6 based upon the redshift 
of arc A and detection in F814W. The other (see Fig.|SJ) 



is Spitzer IRAC-selected with a hkely redshift of z > 6. 
Using this cluster as a gravitational telescope will there- 
fore allow us to study galaxy properties near the epoch 
of reionization. 

The high resolution mass map of this cluster, in addi- 
tion to enabling the high-z galaxy search, will be valu- 
able for several other ongoing science programs. Most 
significantly, we are currently running simulations to re- 
construct the dynamical history of the cluster merger, 
and provide an independent constraint on the dark mat- 
ter self-interaction cross-section. The relative masses of 
the main and the subcluster, and the central mass distri- 
bution, two critical parameters for the simulations, were 
the limiting source of uncertainty prior to this work. 
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Fig. 6. — The effects of different initial conditions (a-c) and exclusion of the strong lensing data (d) on the final reconstruction. Shown 
is the surface mass density k of the weak and strong lensing reconstruction using the two component NIE model (top-left) as initial guess, 
using weak lensing reconstruction from Clowe et al. (2006a) [top-right), and using zero initial model {bottom-left). The weak lensing only 
mass reconstruction is shown in (bottom-right); all for a fiducial source at a redshift of — ♦ oo. The contour levels are linearly spaced with 
Ak = 0.2, starting at k = 0.5. The two cD's and the BCG of the subcluster are denoted with crosses. 
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Fig. 7. — The surface mass density k of the two reconstructions where we {left) removed image system D, and {right) image system I 
from the data. These are two extreme examples where the resulting decreased at most; still the resulting maps are very similar to the 
original reconstruction. The contour levels are the same as in Fig. |^ the two cD's of the main and the BCG of the subcluster are denoted 
with crosses. 




Fig. 8. — The critical curve for a fiducial source at 2 — > oo from the combined strong and weak lensing reconstruction. Circles denote 
the positions of the multiple image systems we used, while the crosses denote the candidate high redshift source identified using Spitzer 
data (note that this system was not included in the reconstruction). North is up and East is left, the field is 4.9' X 3.2'. 
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Fig. 9.— The F435W-F606W-F814W color composite of the 1.1' x 0.9' inner part of the subcluster (left) and of the 0.6' X 0.5' cutout 
around arc B (right) with large circles denoting the positions of the multiple image systems we used. The small circles show the candidate 
counter images as predicted by the weak and strong Icnsing reconstruction (not used in the final reconstruction). Image courtesy of Mischa 
Schirmer. 



